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ABSTRACT 

We investigate the linear growth and vertical structure of the magnetorotational insta- 
bility (MRI) in weakly ionised, stratified accretion discs. The magnetic field is initially 
vertical and dust grains are assumed to have settled towards the midplane, so charges 
are carried by electrons and ions only. Solutions are obtained at representative radial 
locations from the central protostar for different choices of the initial magnetic field 
strength, sources of ionisation, disc structure and configuration of the conductivity 
tensor. 

The MRI is active over a wide range of magnetic field strengths and fluid conditions 
in low conductivity discs. Moreover, no evidence was found of a low-limit field strength 
below which unstable modes do not exist. For the minimum-mass solar nebula model, 
incorporating cosmic ray ionisation, perturbations grow at 1 AU for i? < 8 G. For 
a significant subset of these strengths (200 mG < B < 5 G), the growth rate is 
of order the ideal MHD rate {0.75ft). Hall conductivity modifies the structure and 
growth rate of global unstable modes at 1 AU for all magnetic field strengths that 
support MRI. As a result, at this radius, modes obtained with a full conductivity 
tensor grow faster and are active over a more extended cross-section of the disc, than 
perturbations in the ambipolar diffusion limit. For relatively strong fields (e.g. B > 200 
mG), ambipolar diffusion alters the envelope shapes of the unstable modes, which peak 
at an intermediate height, instead of being mostly flat as modes in the Hall limit are 
in this region of parameter space. Similarly, when cosmic rays are assumed to be 
excluded from the disc by the winds emitted by the magnetically active protostar, 
unstable modes grow at this radius for B < 2 G. For strong fields, perturbations 
exhibit a kink at the height where x-ray ionisation becomes active. 

This study shows that, despite the low magnetic coupling, the magnetic field is 
dynamically important for a large range of fluid conditions and field strengths in 
protostellar discs. An example of such magnetic activity is the generation of MRI 
unstable modes, which are supported at 1 AU for field strengths up to a few gauss. 
Hall diffusion largely determines the structure and growth rate of these perturbations 
for all studied radii. At radii of order 1 AU, in particular, it is crucial to incorporate 
the full conductivity tensor in the analysis of this instability, and more generally in 
studies of the dynamics of astrophysical discs. 
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1 INTRODUCTION 

The magnetorotational instability (MRI; Balbus & Hawley 
1991, 1998; Hawley & Balbus 1991) generates and sustain 
angular momentum transport in differentially rotating astro- 
physical discs. It does so by converting the free energy source 
contributed by differential rotation into turbulent motions 
(e.g. Balbus 2003), which transport angular momentum via 
Maxwell stresses. 

Most MRI models in non-ideal MHD conditions adopt 



either the ambipolar diffusion (Blaes & Balbus 1994, 
MacLow et al. 1995 and Hawley & Stone 1998) or resis- 
tive approximations (Jin 1996, Balbus & Hawley 1998, Pa- 
paloizou & Terquem 1997, Sano & Miyama 1999; Sano, In- 
utsuka & Miyama 1998, Sano et. al. 2000, Fleming, Stone 
& Hawley 2000 and Stone & Fleming 2003). The inclusion 
of Hall conductivity terms is a relatively recent develop- 
ment (Wardle 1999 (W99 hereafter), Balbus & Terquem 
2001, Sano & Stone 2002a,b; 2003, Salmeron & Wardle 2003 
(SW03 hereafter) and Desch 2004). When the Hall effect 
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dominates over ambipolar diffusion, fluid dynamics is de- 
pendent on the alignment of tfio magnetic field with the 
angular velocity vector of the disc (Wardle & Ng 1999), and 
wave modes supported by the fluid are intrinsically modi- 
fied. For example, left and right-circularly polarised Alfven 
waves travel at different speeds and damp at different rates 
in this regime (Pilipp ct al. 1987, Wardle & Ng 1999). Both 
the structure and growth rate of MRI perturbations can be 
substantially modified by Hall conductivity, especially when 
the coupling between ionised and neutral components of the 
fluid is low (W99, SW03). 

In a provions paper wo presented a linear analysis of the 
vertical structure and growth of the MRI in weaJily ionised, 
stratifled accretion discs (SW03). In that study the compo- 
nents of the conductivity tensor were assumed to be constant 
with height. The obtained solutions illustrate the properties 
of the MRI when different conductivity regimes are dom- 
inant over the entire cross-section of the disc. We found 
that when the magnetic coupling is weak, modes computed 
with a non-zero Hall conductivity grow faster and act over a 
more extended cross-section than those obtained using the 
ambipolar diffusion approximation. The height above the 
midplane where the fastest growing modes pealc depends on 
the conductivity regime of the fluid. When ambipolar dif- 
fusion is important, perturbations peak at a higher z than 
when the fluid is in the Hall limit. Furthermore, when the 
coupling is weak, perturbations computed with a full con- 
ductivity tensor peak at different heights depending on the 
orientation of the magnetic field with respect to the angu- 
lar velocity vector of the disc. This is a consequence of the 
dependency of the Hall effect on the sign of f2 ■ B ( W99) . 
Finally, we showed that in the Hall regime perturbations 
can have a very complex structure (high wavenumber) , par- 
ticularly when the magnetic field is weak. In these condi- 
tions, many modes were found to grow, even with a very 
weak magnetic coupling. These results suggest that signifi- 
cant accretion can occur in regions closer to the midplane of 
astrophysical discs, despite the low magnetic coupling, due 
to the large column density of the fluid. This idea contrasts 
with the commonly accepted view that accretion is impor- 
tant primarily in the surface regions, where the coupling be- 
tween ionised and neutral components of the fluid is much 
stronger, but the column density is significantly smaller. 

In a real disc, the components of the conductivity tensor 
vary with height (e.g. Wardle 2003) as a result of changes 
in the abundances of charged particles, fluid density and 
magnetic field strength. This, in turn, is a consequence of 
changes in the ionisation balance within the disc, which re- 
flects the equilibrium between ionisation and recombination 
processes. The former are primarily non-thermal, triggered 
by cosmic rays, x-rays emitted by the central protostar and 
radioactive materials. The later can, in general, occur both 
in the gas phase and in grain surfaces (e.g. Oppenheimer & 
Dalgarno 1974, Spitzer 1978, Umebayashi & Nakano 1980, 
Nishi, Nakano & Umebayashi 1991, Sano et. al. 2000). As a 
result, different conductivity regimes are expected to domi- 
nate at different heights (Wardle 2003). In this paper we re- 
visit the linear growth and structure of MRI perturbations 
using a height-dependent conductivity tensor. We assume 
that the disc is thin and isothermal and adopt a flducial disc 
model based in the minimum-mass solar nebula (Hayashi 
1981, Hayashi, Nakazawa & Nakagawa 1985). Further, we 



assume that ions and electrons are the main charge carriers, 
which is a valid approximation in late stages of accretion, 
after dust grains have settled towards the midplane of the 
disc. As an indication of the timescales for this settling to 
occur, Nakagawa, Nakazawa & Hayashi (1981) show that 
the mass fraction of ~ 1 - 10 pm grains well mixed with the 
gas phase, (i.e. not settled), drops from ~ 10"^ to 10"'* in 
t ~ 2 X 10'' to 1 X 10^ years. Furthermore, although the 
timescale for dust settling all the way to the equator may 
exceed the lifetime of the disc, grains can settle within a few 
pressure scaleheights about the midplane in a much shorter 
timescale (Dullcmond & Dominik 2004). 

This paper is structured as follows: Section 2 sum- 
marises the formulation, including the governing equations, 
fiducial disc model and ionisation balance. Section 3 de- 
scribes the linearisation, parameters of the problem and 
boundary conditions. Results are presented in sections 4 and 
5. We summarise in section 4 the test models used in this 
study and present the ionisation rates as a function of z 
at representative radial locations from the central protostar 
[R = 1, 5 and 10 AU). The importance of different con- 
ductivity regimes at different heights is also described for 
a range of magnetic field strengths. Then, in section 5, the 
structure and growth rate of MRI perturbations at the radii 
of interest arc analysed, including a comparison with results 
obtained using different configurations of the conductivity 
tensor, sources of ionisation and disc structure. We find that 
the MRI is active over a wide range of fluid conditions and 
magnetic fleld strengths. For example, for the flducial model 
at i? = 1 AU, and including cosmic ray ionisation, unstable 
modes are found for B < 10 G. When 200 mG < B < 5 
G, these perturbations grow at about the ideal-MHD rate 
~ 0.75 times the dynamical (keplerian) frequency of the disc. 
Results are discussed in section 6 and the paper's formula- 
tion and key findings are summarised in section 7. 



2 FORMULATION 

2.1 Governing Equations 

Following W99 and references therein, we formulate the con- 
servation equations about a local keplerian frame corotating 
with the disc at the keplerian frequency Q = ^ GM/r^. 
Time derivatives in this frame, d/dt, correspond to d/dt + 
0.8/ d4> in the standard laboratory system {r,(f>,z) anchored 
at the central mass M. The fluid velocity is expressed as a 
departure from keplerian motion v = V — vk, where V is 
the velocity in the laboratory frame and vk = ^jGMjr (f) 
is the keplerian velocity at the radius r. The fluid is as- 
sumed to be weakly ionised, so the effect of ionisation and 
recombination processes on the neutral g £IS, £LS well as the 
ionised species' inertia and thermal pressure, are negligible. 
Under this approximation, separate equations of motion for 
the ionised species are not required and their effect on the 
neutrals is treated via a conductivity tensor (W99 and ref- 
erences therein), which is a function of location {r,z). 
The governing equations are the continuity equation, 

g+V-(pt;) = 0, (1) 
the equation of motion. 
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and the induction equation, 
dB 



dt 



= V X (v X B) - cV X - ^nBr^ . 



(2) 



(3) 



In the equation of motion (2), $ is the gravitational potential 
of the central object, given by 

GM 



(r2 + 22)2 



(4) 



and v%/r is the centripetal term generated by exact kep- 
lerian motion. Similarly, 20v^f and ^Q,Vr<^ are the coriohs 
terms associated with the use of a local keplerian frame, 
Cs = ^/P/p is the isothermal sound speed, Q = vx/r is 
the keplerian frequency and c is the speed of light. Other 
symbols have their usual meanings. 

In the induction equation (3), E' is the electric field in 
the frame comoving with the neutrals and the term ^QBr4> 
accounts for the generation of toroidal field from the radial 
component due to the differential rotation of the disc. 

Additionally, the magnetic field must satisfy the con- 
straint: 

V • B = , (5) 
and the current density must satisfy Ampere's law. 



J= -^V X B 

and Ohm's law, 

J = c ' E' = <^||-^fi -\- (T\B Ej_ + o'2Ej_ . 



(6) 



(7) 



Note that we have introduced the conductivity tensor cr 
in equation (7). We refer the reader to W99 and references 
therein for derivations and additional details of this formu- 
lation. Assuming that the only charged species are ions and 
electrons, and that charge neutrality is satisfied {m = Ue), 
the components of the conductivity tensor can be expressed 
as (SW03), the conductivity parallel to the magnetic field, 

ail = 2^ (ft - (8) 

the Hall conductivity, 

B + + ' 
and the Pedersen conductivity, 

B (1-f + ■ 
Prom (9) and (10), we find: 
eerie (A - A) 



(9) 



(10) 



B [(l+/32)(l+;32)]l/2 



(11) 



where crx — Vof + of is the total conductivity perpendic- 
ular to the magnetic field. 
In equations (8) to (11), 



ZjeB 1 



(12) 



is the Hall parameter, given by the ratio of the gyrofrequency 
and the collision frequency of charged species j with the 
neutrals. It represents the relative importance of the Lorentz 
and drag terms in the charged species' motion. In equation 
(12), 



li = 



< av >j 
ruj +m 



(13) 



where m is the mean mass of the neutral particles and 

< (TV >j is the rate coefficient of momentum exchange by 
collisions with the neutrals. The ion-neutral momentum rate 
coefficient is given by. 



< av >i= 1.6 X 10 



) 3 -1 

cm s 



(14) 



an expression that ignores differences in the values of elastic 
cross-sections of H, H2 and He (Draine, Roberge & Dalgarno 
1983). Similarly, the rate coefficient of momentum exchange 
of electrons with the neutrals is (Draine et al. 1983): 



<a.>..lxlO-ri^V'%m-. 



(15) 



and the mean ion mass rrii = SOttih- 

The relative magnitudes of the components of the con- 
ductivity tensor differentiate three conductivity regimes: 
When (T|| ^ (72 ^ |a"i| for most charged species, ambipo- 
lar diffusion dominates and the magnetic field is effectively 
frozen into the ionised components of the fluid. Electron-ion 
drift is small compared with ion-neutral drift and the ionised 
species act as a single fluid. MRI studies in this regime in- 
clude Blaes & Balbus (1994), MacLow et al. (1995) and Haw- 
ley & Stone (1998). Conversely, when k, 02 ^ Icil, the 
conductivity is a scalar, giving rise to the familiar ohmic dif- 
fusion, and the magnetic field is no longer frozen into any 
fluid component. Examples of studies of the MRI in this 
regime are Jin (1996), Balbus & Hawley (1998), Papaloizou 
& Terquem (1997), Sano & Miyama (1999), Sano et. al. 
(1998, 2000, 2004), Fleming, Stone & Hawley (2000) and 
Stone & Fleming, (2003). Ambipolar diffusion dominates in 
low density regions, where magnetic stresses axe more im- 
portant than collisions with the neutrals and the ionised 
species are mainly tied to the magnetic field rather than to 
the neutral particles. On the contrary, in relatively high den- 
sity environments, the ionised species are primarily linked to 
the neutrals via collisions and Ohmic diffusion dominates. 
There is, however, an intermediate density range charac- 
terised by a varying degree of coupling amongst charged 
species. In these circumstances, there is a component of the 
conductivity tensor that is perpendiculax both to the elec- 
tric and magnetic fields. It is this term that gives rise to Hall 
currents. It has been shown that this regime can be impor- 
tant in the weakly ionised environment of accretion discs. 
For example, using an MRN grain model (Mathis, Rumpl 
& Nordsieck 1977) with a power law distribution of grain 
sizes between 50 and 2500 Angstrom, Hall conductivity is 
important for lO'^ < nn < lO" cm"^ (Wardle & Ng 1999). 
MRI studies including Hall diffusion have been conducted 
by W99, Balbus & Terquem (2001), Sano & Stone (2002a,b; 
2003), SW03 and Desch (2004). 

The ionisation balance within the disc (section 2.3) de- 
termines the abundances of charged species (ions and elec- 
trons). These, in turn, determine the configuration of the 
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conductivity tensor. In protostellar discs, outside of the cen- 
tral 0.1 AU, iouisation sources arc non-thermal (Hayashi 
1981) and the ionisation fraction is not enough to produce 
good magnetic coupling over the entire cross-section of the 
disc. In these conditions, the region around the midplane is 
likely to be a magnetically "dead zone" (Gammie 1996). 



2.2 Disc Model 

Our model incorporates the vertical structure of the disc, 
but neglects fluid gradients in the radial direction. This is an 
appropriate approximation, as astrophysical accretion discs 
arc generally thin and changes in the radial direction occur 
on a much bigger length scale than those in the vertical di- 
rection. Including the vertical structure means that pertur- 
bations of spatial dimensions comparable to the scale height 
of the disc, which are associated with either a strong mag- 
netic field {va ~ Cs) or low conductivity, can bo explored. 

We adopted, as our fiducial model, a disc based in the 
minimum-mass solar nebula (Hayashi 1981, Hayashi et. al. 
1985). In this model, the surface density S(r), sound speed 
Cs(r), midplane density Po(r), scaleheight H{r) and temper- 
ature T(r) are: 



^0 



Po(r)=p„(^)" 

^W = ^°(iIZ7) ' 



2.34 



and, 



1/4 



In the previous expressions. 
So = 1.7 X lO^gcm"^, 
Cao = 9.9 X 10* cm s"^ , 
Po = 1.4 X 10"®g cm"^ , 
Ho = 5.0x 10"cm, 
and, 

To = 280K. 



(16) 
(17) 

(18) 
(19) 

(20) 

(21) 
(22) 
(23) 
(24) 

(25) 



In equations (16) to (20), Af* and are the stellar mass 
and luminosity, respectively, and /i is the mean molecular 
mass of the gas. This model describes the minimum mass dis- 
tribution of the solax nebula, estimated assuming an efficient 
planet formation with no significant migration. With these 
assumptions, the current mass distribution and composition 
of the planets is a good indication of that of dust in the orig- 
inal nebula. This model has been used extensively, but there 
are theoretical grounds to expect that a typical protostellar 
disc may be more massive, with a different large scale struc- 
ture, shaped ultimately by the action of MHD turbulence 
(e.g. Balbus & Papaloizou 1999). It can be shown that the 



disc surface density may be roughly up to an order of mag- 
nitude higher than the one specified by the minimum-mass 
solar nebula model before self gravity becomes important. 
To account for the possibility of discs being more massive 
than the minimum-mass solar nebula model, we also studied 
the properties of MRI unstable modes using a disc structure 
with an increased surface density = lOSo and mass den- 
sity Po = lOpc,. The Toomre parameter (Toomre 1964) Q ~ 6 
in this case, so the assumption of a non self-gravitating disc 
is still valid. For simplicity, we assume that the tempera- 
ture To is unchanged, so Ho and Cao are the same as in the 
fiducial model. 

As the disc is gravitationally stable, the gravitational 
force in the vertical direction comes from the central proto- 
star. Under these conditions, the balance between the ver- 
tical component of the central gravitational force and the 
pressure gradient within the disc determines its equilibrium 
structure. The vertical density distribution in hydrostatic 
equilibrium is given by 



po{r) 



— exp 



2/f2(r) 



(26) 



Assuming a neutral gas composed of molecular hydrogen 
and helium such that uhb = 0-2nH2 > the neutral gas mass 
and number densities are, respectively. 



p = Hmmi 
and 



lAmHTlH 



n = 1.2nH2 I 
which gives, 

p(r, z) 



nH{r,z) 



lAruH 



(27) 



(28) 



(29) 



and ^ = p/n = 2.34. For simplicity, we take Li,/LQ = 
Mi,/MQ = 2.34/yti = 1 in all our models. 

2.3 Ionisation balance 

The ionisation balance within the disc is given by the equi- 
librium between ionisation and recombination processes. In 
general, recombination can take place both in the gas phase 
and on grain surfaces (e.g. Oppcnhcimer & Dalgarno 1974, 
Spitzer 1978, Umebayashi & Nakano 1980, Nishi, Nakano 
& Umebayashi 1991, Sano et. al. 2000), but here we have 
assumod that grains have settled, so wo arc including only 
gas-phase recombination rates. Except in the innermost sec- 
tions of the disc (i? < 0.1 AU), where thermal efi^ects are im- 
portant (Hayashi 1981), ionisation processes in protostellar 
discs are mainly non-thermal. Ionising agents are typically 
cosmic rays, x-rays emitted by the magnetically active proto- 
star, and the decay of radioactive materials present within 
the disc. Some authors (e.g. Promang, Terquem & Balbus 
2002) have argued that the low energy particles important 
for cosmic-ray ionisation are likely to be excluded by the 
protostar's winds. Given the uncertainties involved, and the 
expectation that cosmic rays (if present), may be more im- 
portant than x-rays near the midplane (especially for i? w 1 
AU, where the surface density is larger than the attenua- 
tion length of x-rays) , we explore both options in this study. 
The treatment of ionisation and recombination processes is 
detailed below. 
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2.3.1 Cosmic ray ionisation 

The cosmic ray ionisation rate, CcR{r,z), is given by, 



CcR{r, z) = ^ 



exp - 



^CR 



exp 



S(r)-S(r,z) 

Af7R 



(30) 



where (^cr = W~ is the ionisation rate due to cosmic 

rays in the interstellar medium and Xcr = 96 gr cm~^ is 
the attenuation length, a measure of how deep cosmic rays 
can penetrate the disc (Umebayashi & Nakano 1981). When 
the surface density is larger than ~ 2Xcr, most cosmic rays 
do not reach the midplane. Also, 



p oc 

E(r, z)= p{r, z)dz , 

•/ Z 



(31) 



is the vertical column density from the location of interest 
to infinity. The two terms in brackets in equation (30) mea- 
sure the ionisation rate at the position (r, z) by cosmic rays 
penetrating the disc from above, and below, respectively. 

2.3.2 X-ray ionisation 

There is strong evidence for an enhanced magnetic activ- 
ity in young stellar objects (i.e. sec review by Glassgoki, 
Feigelson & Montmerlc 2000). Typical soft x-ray luminosi- 
ties (0.2 - 2 kev) of these objects arc in the range of 10'^*^ - 
lO^^ergs"^, about 10^ - 10^ times more energetic than solar 
levels (Glassgold et. al. 2000). Models of the penetration of 
stellar x-rays into a protostellar disc by Igea & Glassgold 
(1999), show that even discounting the low energy photons 
that are attenuated by stellar wiuds, the ionisation rate due 
to hard x-rays close to the central object is many orders of 
magnitude larger than that of cosmic rays, especially above 
z/H ~ 2. These authors investigated the ionisation rate by 
x-rays from a central protostar modelled as a x-ray source 
of total luminosity Lx ~ lO^^ergs^^ and temperature kT^ 
in the range of 3 - 8 kev. The transport of x-rays through 
the disc was followed using a Monte Carlo procedure which 
included both absorption and scattering by disc material. 
The incorporation of scattering is important, as x-rays can 
be scattered not only out of the disc, but also towards the 
midplane, enhancing the ionisation level deeper within the 
disc. Results indicate that, at each radius, the x-ray ioniser 
tion rate is a function of the vertical column density into the 
disc N±{r,z), irrespective of its structural details. 

To calculate the x-ray ionisation rate Cx {r, z) in the up- 
per half of the disc we added the contribution of x-rays arriv- 
ing from both sides. We used the values of Cx(r, z) ( s~^) as a 
function of the vertical column density, N±_ (cm~^), plotted 
in Fig. 3 of Igea & Glassgold (1999), for ii = 1, 5 and 10 AU 
and kTx = 5 Kev. The vertical column density, appropriate 
for x-rays arriving from the top is 



iV±(r,2) 



-f 



(32) 



with uh given by equation (29). Similarly, the ionisation 
contributed by x-rays arriving at (r, 2; > 0) from the other 
hemisphere of the disc is obtained substituting N± above by 
2A''x(r, 0) — N±{r,z), although this contribution is usually 
negligible. 



2.3.3 Radioactivity 

The ionisation rate C« contributed by the decay of radioac- 
tive materials (primarily *^K) have also been considered in 
previous work (i.e. Consolmagno & Jokipii 1978, Sano et 
al. 2000). This rate can be calculated as (Consolmagno & 
Jokipii 1974): 



Cr' 



-yUrE 



(33) 



where 7 ( s~^) is the decay rate of the radioactive isotopes, 
Ur their number density relative to hydrogen, and E (eV) the 
energy of the produced radiation. Similarly, Ae = 37 eV, is 
the average energy for the production of an ion pair in H2 gas 
(Consolmagno & Jokippi 1974 and references therein, ShuU 
& Van Steenburg 1985, Voit 1991). In the present study we 
have adopted, 



Cfi = 6.9 X lO^'^^lJa (1 - 52) fg] s- 



(34) 



(Umebayashi & Nakano 1981, 1990). Here, 82 « 0.02 is the 
fraction of heavy metal atoms in the gas phase, estimated 
via measurements of interstellar absorption lines in diffuse 
clouds (Morton 1974). On the other hand, /g is a parameter 
that takes into account the degree of sedimentation of dust 
grains in protostellar discs with respect to interstellar values 
(Hayashi 1981, Sano et. al. 2000). 

Although the ionisation effect contributed by this agent 
is very small compared with that of x-rays and cosmic rays, 
we included it because it may well be the only mechanism ac- 
tive in regions close to the midplane (particularly for i? < 5 
AU) in the scenario where cosmic rays are excluded from the 
disc. In this case, it is important to explore the sensitivity 
of MRI perturbations to changes in the level of depletion of 
dust grains (/g). 

We have assumed here that the ionised component of 
the fluid is made of ions and electrons only, a case that corre- 
sponds to late evolutionary stages of protostellar discs, after 
dust grains have settled towards the midplane. Accordingly, 
we assume /g = in all our models, except when especifi- 
cally exploring the role of radioactivity in a disc where cos- 
mic rays are assumed to be excluded from it. 



2.3.4 Recombination rate 

Gas-phase recombination occurs through dissociative recom- 
bination of electrons with molecular ions and, at a slower 
rate, via radiative recombination with metal ions. It has 
been pointed out by previous authors that the ionisation 
balance may be especially sensitive to the presence of metal 
atoms within the disc. For some disc configurations, a num- 
ber density of metals as small as lO"'^ times the cosmic 
abundance may be enough to make turbulent the whole 
cross-section of the disc, eliminating the central magnetic 
dead zone (Fromang et. al. 2002). This is so because metal 
atoms generally take the charges of molecular ions, but re- 
combine with electrons (via radiative processes) at a much 
more slower rate. 

If and the number densities of molecular 

and metal ions, respectively, the rate equations for rie and 
n+ can be expressed as (e.g. Fromang et. al. 2002): 



dUe ^ 



(35) 
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and 



3 METHODOLOGY 



dnt, 
dt 



+ 

m ? 



(36) 



where C is the total ionisation rate, calculated as suriirnarised 
in the previous sections, and nn is the hydrogen number 
density from equation (29). 

In the previous expressions, /3 is the dissociative recom- 
bination rate coefficient for molecular ions, /3r is the radia- 
tive recombination rate coefficient for metal ions and /3t is 
the rate coefficient of charge transfer from molecular ions 
to metal atoms. If all metals are locked into dust grains, 
which have in turn sedimontcd towards the midplane of the 
disc, equations (35) and (36), together with charge neutral- 



ity, rie = r 
coefficients. 



= 3 X 10" 



+ rQn, and appropriate values for the rate 



-1/2 3 - 

' cm s 



/3 = 3x IQ-^T-'/^cm^-V 



/3t = 3 X 10"^cm^s"\ 



(37) 
(38) 
(39) 



(see Promang et al 2002 and references therein) , lead to 



(,nH 



(40) 



However, as these authors point out also, dust grains not 
only absorb metal atoms, but also release them as a result 
of the action of x-rays. Because of this effect, the abundance 
of metal atoms in a disc could be quite insensitive to the 
spatial distribution of dust grains. This, in turn, means that 
dust settling does not necessarily lead to a severe reduction 
in the number density of metal atoms in the gas phase. If 
these dominate, the corresponding rie is 



/3r 



(41) 



The transition from one regime to the other may be shown 
to occur when (Fromang et. al. 2002), 



Xm 



10"^T" 



1/2 



(42) 



It is clear then, that the evolutionary stage and 
activity- of the disc are important factors in the ionisation 
balance, as they influence the degree of sedimentation of 
dust grains. As the disc evolves, dust grains tend to occupy 
a thin layer around the midplane, becoming removed -in any 
dynamical sense- from the gas at higher vertical locations. 
This causes the ionisation fraction of the gas to increase, 
by eliminating recombination pathways on dust surfaces. In 
the present work, we have used equation (41) to calculate 
the electron (and ion) number density. The minimum val- 
ues of Xm for this approximation to be valid (equation (42)) 
have also been computed and compared with an estimate of 
metal abundances in the gas phase (see section 4.2). Results 
indicate that, for the range of parameters adopted here, the 
abundance of metal atoms is such that radiative recombi- 
nation is indeed dominant, so the use of equation (41) to 
calculate the electron fraction is justified. 



3.1 Linearisation 

Pull details of the methodology are described in SW03. For 

the sake of clarity, and completeness, we summarise here the 
most important stops and point at some differences with 
the previous formulation. The system of equations (1) to 
(3), (6) and (7) was linearised about an initial steady state 
where fluid motion is exactly keplerian and the magnetic 
field is vertical, so J = v = E' = Q and B — Bz. In the 
initial state both E' and J vanish, so the changes in the 
conductivity tensor due to the perturbations do not appear 
in the linearised equations. As a result, it is not necessary 
to explore how the perturbations affect the conductivity of 
the fluid and only the values in the initial steady state are 
required. 

We assume that the wavevector of the perturbations 
is perpendicular to the plane of the disc {k = kz)- These 
perturbations, initiated from a vertical magnetic field, are 
the fastest growing modes when the fluid is in either the 
Hall or resistive conductivity regimes, as in these cases mag- 
netic pressure suppresses displacements with a non-zero rar- 
dial wavenumber (Balbus & Hawley 1991, Sano & Miyama 
1999). However, as pointed out by Kunz & Balbus (2004) 
through their analysis of local axisymmetric MRI, this is not 
necessarily the case for perturbations in the ambipolar dif- 
fusion limit. Their formulation includes radial and vertical 
wavcnumbcrs as well as vertical and azimuthal field compo- 
nents. Results indicate that in this regime, the fastest grow- 
ing modes may exhibit both radial and vertical wavenum- 
bers. Interestingly, in this configuration unstable modes are 
present even for outwardly increasing angular velocity pro- 
files, contrary to the case in ideal-MHD. 

Taking perturbations of the form q = + Sq{z)e^'^* 
about the initial state, linearising and neglecting terms of 
order H/r or smaller, we find that the final linear system of 
equations that describes the MHD perturbations within the 
disc is, 



iujpSvr — 2pQ,5vs 5Ja, = , 

c 

iujpSv^, + ^p^lSVr -\ SJr = , 

c 

. ^ dSE^ 
lUJOBr — c — --^ = , 
dz 



iu>5B^ + c^^ + ln5Br = 0, 

dz ^ 



6Jr = — 
SJ4, 



dSEr 

dz 
c dSB^ 
4n dz 
c dSBr 



47r dz ' 



5Jr = cr25E'^ — aiSE'j, , 



SJ^ = aiSE'r + 025E'j^ 



4> 1 



(43) 
(44) 
(45) 
(46) 
(47) 
(48) 

(49) 

(50) 



where 5E^ and SEr are the perturbations of the electric field 
in the laboratory frame. 



5E^ = 6E'^ + —5vr , and 



(51) 
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SEr =SE'- — 



5vi 



(52) 



Note that , the component of the conductivity tensor 
parallel to the magnetic field, does not appear in the lin- 
earised equations. This, in turn, implies that in the linear 
phase of the MRI and under the adopted approximations, 
the ambipolar diffusion and resistive conductivity regimes 
behave identically (see section 2.1). 

We express equations (43) to (52) in dimensionless form 
by normalising the variables as follows: 



z 
H 



Pir,z) 
Po(r) 



SB* 



5B 

Bo 



5v* = 



5v 



SE* = 



cSE 

CsBo 



SE'* = 



cSE' 

CsBq 



A2 = h vx 



1 



and 



^^ = ^ + 2^1 + i.2- 
In the above expressions, 

Bo 



is the local Alfven speed in the disc, and 



_ Wc _ 1 -Bp ax 

^~ n ~ n pc^ 



(61) 
(62) 

(63) 
(64) 



controls the local coupling between the magnetic field and 

the disc (W99, SW03, see also section 3.2). Non-ideal ef- 
fects strongly modify wavemodes at, or above, the critical 
frequency uic- 



SJ* = 



cSJ 



CsBoa± a± 

Here, subscript 'o' denotes variables at the midplanc 
of the disc. Note that wc have used the local a± instead 
of (T±g, as in SW03, to normalise cr and SJ. This is more 
useful when dealing with a height-dependent conductivity. 
In the following dimensionless system, we have dropped the 
asterisks and expressed the final equations in matrix form: 



dz 



f Br \ / 







ClAl C1A2 \ / Br \ 



-C1A2 C1A3 
-I/O 

y 



B^ 



(53) 



/ A2 -A3 \ 
5 J = C2 \5E 
\ Al A2 



Sv = x 



SJ 




where 



C2 = 



1 + X 



5 cTi 0-2 

,2 I n—+^^—+X 



1 + \2cr± 
1 



a± 1 + 1/^ 



(54) 
(55) 

(56) 

(57) 
(58) 

(59) 
(60) 



3.2 Parameters 

The following parameters control the dynamics and evolu- 
tion of the fluid: 

(i) va/cs, the ratio of the local Alfven speed to the 
isothermal sound speed of the gas. It is a measure of the 
strength of the magnetic field. In ideal-MHD conditions, un- 
stable modes grow when the magnetic field is subthermal 
{va/cs < 1). Under this approximation, when va ~ Cs the 
minimum wavelength of the instability is of the order of the 
scale height of the disc and the growth rate decreases rapidly. 
This is also the case under the assumption that ambipolar 
or ohmic diffusion dominates over the entire cross-section 
of the disc (SW03). However, when the fluid is in the Hall 
regime, with the magnetic field counteraligned with the an- 
gular velocity vector of the disc (f2 ■ B < 0), MRI unstable 
modes may exist for stronger fields, up to va/cs < 3 (SW03). 

(ii) Xi a parameter that characterises the strength of the 
local coupling between the magnetic field and the disc (equa- 
tion (64)). It is given by the ratio of the critical frequency 
above which fiux-freezing conditions break down and the dy- 
namical (keplerian) frequency of the disc. If x = Wc/f2 < 1, 
the disc is poorly coupled to the disc at the frequencies of in- 
terest for dynamical analysis, which are also the interesting 
frequencies for the study of the MRI. 

(iii) (Ji/a2, the ratio of the conductivity terms perpendic- 
ular to the magnetic field. It is an indication of the conduc- 
tivity regime of the fluid, as discussed in section 2.1 . 

We calculated the values of the three parameters above 
at different locations (r, z) within the disc using equations 
(9), (10), (17), (63) and (64). The strength of the magnetic 
fleld is a free parameter of the problem. 



3.3 Boundary Conditions 

To solve equations (53) to (56) it is necessary first to inte- 
grate the system of ordinary differential equations (ODE) 
in (53). This can be treated as a two-point boundary value 
problem for coupled ODE. Five boundary conditions must 
be formulated, prescribed either at the midplane or at the 
surface of the disc. 



© 2004 RAS, MNRAS 000, 1-26 



8 R. Salmeron and M. Wardle 



At the midplane: We chose to assign odd symmetry to 

5 Br and SB^, so they vanish at z = 0. This gives us two 
boundary conditions at the midplane. Also, as the equa- 
tions are linear, their overall scaling is arbitrary, and a third 
boundary condition can be obtained by setting one of the 
fluid variables to any convenient value. To that effect, we 
assigned a value of 1 to (5_E^. The three boundary conditions 
applied at the midplane are, then: 

SBr = 5B^ = 0, and 

SE'r = 1 . 

At the surface: x is inversely proportional to the den- 
sity (equation 64), so if the conductivity is assumed to be 
constant with height, it increases monotonically with z. This 
was the case in SW03, where we used this argument to pro- 
pose that at sufficiently high z above the midplane, ideal 
MHD conditions held. This, in turn, implies that the wave- 
lengths of magnetic field perturbations, given the adopted 
dependency of p with z (equation (26)), must tend to infin- 
ity when z ^ oo. As a result, the amplitude of such modes 
should vanish at infinity as well. Here, however, the con- 
ductivities vary with so this argument needs revisiting. 
In section 4.3 we analyse the dependency of x with height 
for different radial positions and a range of magnetic field 
strengths. Here, we just highlight that at the surface of the 
disc {z/H = 6, see below), the magnetic coupling is still 
strong. For example, for the fiducial model at = 1 AU, 
X ~ 20. It decreases to ~ 4.5 for i? = 5 AU and to ~ 2.5 
for _R = 10 AU. This is still above the limit (~ 1) for strong 
coupling (W99). 

When X > 10 the growth rate and characteristic wavenum- 
ber of local unstable modes differ little from the ideal limit 
(W99), so the same line of reasoning of SW03 can be used 
to argue that in this case SBr and 5B^ should be zero at the 
boundary as well. Although x at the surface when -R > 5 
AU is below 10, it is considered that the magnetic coupling 
there is still sufficiently strong for the ideal-MHD approxi- 
mation to be essentially valid as well. This is also confirmed 
by the way in which magnetic field perturbations tend to 
zero at the boundary at these radii (see results in sections 4 
and 5). Consistently with this, both 5E' and 5 J vanish at 
the surface as well. This gives us the remaining two bound- 
ary conditions required to integrate the system of equations 
(53). 

We chose to locate the boundary at z/H = 6, after confirm- 
ing that increasing this height does not significantly affect 

either the structure or the growth rate of unstable modes. 
Summarising, the boundary conditions adopted at the sur- 
face are: 

SBr = 5B^ = 

at z/H = 6. 

This system of equations is solved by 'shooting' from 
the midplane to the surface of the disc while simultaneously 
adjusting the growth rate v and SE'^. The adjustment is 
done via a multidimensional, globally convergent Newton- 
Raphson method, until the solution converges. 



4 DISC CONDUCTIVITY 
4.1 Test Models 

Our fiducial model is based in the minimum-mass solar neb- 
ula disc (Hayashi 1981, Nakazawa & Nakagawa 1985; see 
section 2.2). The structure and growth rate of MRI unsta- 
ble modes are calculated at representative radial positions 
{R = 1, 5 and 10 AU) from the central protostar. Two sce- 
narios are explored: Cosmic rays either penetrate the disc, 
although attenuated as appropriate (as given by eq. (30)), or 
they are excluded from it by the winds produced by the cen- 
tral object. Unless stated otherwise, results presented using 
this model include cosmic ray ionisation. We also consider 
a disc model with an increased mass and surface density, 
as described in section 2.2. For simplicity, other disc pa- 
rameters remain unchanged and in this case, all calculations 
incorporate the ionisation rate provided by cosmic rays. 

We also compare solutions obtained using different con- 
figurations of the conductivity tensor. By comparing the 
structure and growth rate of unstable modes found using 
the commonly adopted ambipolar diffusion -or resistive- 
(cTi = 0) approximations with those obtained with a full 
conductivity tensor {aiBz > 0, a2 7^ 0), we can appreciate 
how, and in which regions of parameter space, Hall conduc- 
tivity alters the properties of the instability. This is explored 
for the fiducial model at 1 AU only. 

The properties of the MRI in the Hall limit are depen- 
dent on the alignment of the magnetic field and angular 
velocity vectors of the disc ( W99) . The case when these vec- 
tors are parallel (antiparallel) is characterised by aiB^ > 
{(JiBz < 0). In the Hall {aiBz < 0) limit, our code fails to 
converge whenever the combination of parameters is such 
that X < 2 anywhere in the domain of integration. This 
is not surprising as, in this regime, when 0.5 < x < 2 all 
wavenumbers are unstable (W99). As a result, we explored 
the effect of the alignment of B and f2 by compaxing solu- 
tions obtained with a full conductivity tensor but incorpo- 
rating o\Bz terms of opposite sign. Even with this approach, 
solutions including a. aiB^ < term could not be computed 
at all for i? = 1 AU, given that the midplane coupling at 
this radius is extremely low (~ 10"^° to 10~^ for B between 
1 mG and 1 G). Therefore, we present this analysis at 5 and 
10 AU only. Solutions incorporating a negative Hall conduc- 
tivity could be obtained for B > 8 mG at 5 AU and B > 1.5 
mG at 10 AU, where the magnetic coupling at the midplane 
is > 2. 

Unless stated otherwise, results discussed in the follow- 
ing sections incorporate a,a\Bz > Hall term contribution. 



4.2 Ionisation Rates 

Fig. 1 shows the ionisation rate contributed by cosmic rays, 
x-rays and radioactive materials as a function of height for 
the fiducial model. Results are presented at the radial loca- 
tions of interest: 1, 5 and 10 AU. Note that x-rays penetrate 
up to z/H ~ 1.7 and z/H ~ 0.3 for 7? = 1 and 5 AU, re- 
spectively. They reach the midplane for i? = 10 AU. The 
x-ray ionisation rate is heavily attenuated below ~ 1.5 to 2 
scaleheights from the surface, depending on the radius. 

At 1 and 5 AU, the regions where x-rays are not able 
to penetrate are ionised only by the action of cosmic rays 
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(if present) and radioactive decay. Moreover, cosmic rays 
constitute the most efficient ionising agent for z/H < 2.2 
at all studied radii. As a result, if they are excluded from 
the disc by the protostar's wiuds (i.e. Fromang et al. 2002), 
the magnetic coupling of the gas in these sections of the 
disc is expected to be severely reduced. In order to explore 
the properties of MRI perturbations under this assumption, 
we compare results computed including and excluding this 
source of ionisation for the fiducial model (section 5). 

Igea & Glassgold (1999) have pointed out that in proto- 
stellar discs, the x-ray ionisation rate can be 10^ - 10® times 
that of interstellar cosmic rays. Their calculations show that 
x-rays dominate over cosmic rays until the vertical column 
density, N±, is 10^®, 5 x 10^^ and 2 x 10^^ cm"^ for R=l, 
5 and 10 AU, respectively. Once N± exceeds ^ lO'^'cm"^, 
x-ray ionisation dies out because there are no remaining inci- 
dent photons available. The ionisation rates adopted in this 
study (Fig. 1) are in agreement with these findings. 

We have used eq. (41) to calculate the electron num- 
ber density (section 2.3.4). This expression is valid when 
metal atoms dominate and electrons rccombine primarily 
via radiative processes. Here, we justify this choice. Fig. 2 
compares the minimum abundance of metals, (eq. (42)), for 
this approximation to be valid with an estimate of total 
metal abundances in the gas phase. Minimum abundances 
(relative to hydrogen) are shown for the radial locations of 
interest for the fiducial model. The total gas-phase metal 
abundance is taken as 8.4 x 10~®52 (Umebayashi & Nakano 
1990). From Fig. 2 it is clear that the abundance of gas- 
phase metal atoms is expected to exceed the mininmm for 
radiative recombination processes to be dominant at all radii 
of interest here. As a result, we conclude that equation (41) 
is a valid approximation to calculate the electron number 
density in this study. 



4.3 Magnetic coupling and conductivity regimes 
as a function of height 

It is interesting to explore which conductivity terms (section 
2.1) are dominant at different heights at the representative 
radial locations we are considering. This information will be 
used in the analysis of the structure and growth of MRI 
unstable modes in the next section. We recall that Hall dif- 
fusion is locally dominant when (SW03), 



■10 



X < Xcrit 



|ai| 



(65) 



When X is higher than this, but still < 10, ambipolar 
diffusion dominates. In regions where x is stronger, ideal- 
MHD conditions hold ( W99) . We remind the reader than in 

the present formulation, the ambipolar and Ohmic conduc- 
tivity regimes behave identically. This is the case because 
fj|| (the conductivity component that distinguishes between 
them), does not appear in the final, linearised system of 
equations (see section 3.1). This should be kept in mind 
when analysing the results, as even though the case where 
ai = has been referred to as the 'ambipolar diffusion' limit 
(cry S> (72), it is also consistent with the 'resistive' regime 

(CT|| ~ (72). 

We are interested in comparing ^ and 'Xcrit as a func- 
tion of height for different disc models, radial locations and 
choices of the magnetic field strength. An example of the 
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Figure 1. Ionisation rate ( s~^) contributed by cosmic rays (curve 
labeled 'cr'), x-rays (xr) and radioactive materials (rad) as a func- 
tion of fieiglit, for tfie fiducial model (minimum mass solar nebula 
disc). Results are sliown for R = 1, 5 and 10 AU. Note tfiat at 
1 AU, cosmic ray ionisation is attenuated witfi respect to tfie in- 
terstellar rate for z/H < 2. Also, x-rays are excluded from the 
disc for z/H < 1.7 (at~l AU) and z/H < 0.3 (at 5 AU). For 
R = 10 AU, they reach the midplane, although their ionisation 
rate is heavily attenuated. Cosmic rays (if present) are the most 
efficient ionising agent for z/i? < 2.2 in all cases. 



typical dependency of these variables with z is shown in 
Fig. 3 for the fiducial model at il = 1 AU for B = 10 
mG. Note that the conductivity term parallel to the field 
increases monotonically with z and is independent of B (see 
eqs. (8) and (12)). In this case, ai is nearly two orders of 
magnitude smaller than fT2 at the midplane. As a result, near 
the midplane (7|| ~ (T2 » cri, and the fluid is in the resis- 
tive regime (section 2.1). This is consistent with findings by 
Wardle (in preparation), that Ohmic conductivity is impor- 
tant at this radius for relatively weak fields {B < 100 mG) 
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O 




Figure 2. Comparison of tlic minimum abundance of metal 
atoms in the gas pliase, xm (n/nu), for tlic radiative recombi- 
nation rate of metal ions (fir) to be the dominant recombination 
mechanism (Fromang et. al. 2002), and an estimate of the total 
gas-phase metal abundance, as a function of height. Minimum 
abundances are shown for the fiducial model at 1 AU (solid line), 

5 AU (dashed line) and 10 AU (dotted line). The estimated total 
gas-phase metal abundance (dot-dashed line) is from Umebayashi 

6 Nakano (1990). Tiiis abundance exceeds the minimum required 
for all radii and vertical locations of interest in this study. 



when the density is sufficiently high (e.g. uh ^ 10^^ cm~^ 
for B = 1 mG). Both components of the conductivity tensor 
perpendicular to the field increase witfi ficight initially, re- 
flecting the enhanced ionisation fraction at higher z. There 
is then a central section of the disc where ai ^ 02 (the fluid 
is in the Hall regime) and finally, both terms drop abruptly 
closer to the surface as a consequence of the fall in fluid den- 
sity. Note that the Hall conductivity term decreases much 
more sharply than the ambipolar diffusion component. As a 
result, the ambipolar diffusion term is typically several or- 
ders of magnitude greater than the Hall term in the surface 
regions of the disc. To reiterate this concept, we repeat here 
that even in regions where |cri | <^ (72, the growth and struc- 
ture of the MRI in the linear regime are still dominated by 
Hall diffusion if x < Wi\/a^ (SW03). 

The previous results are in agreement with the notion 
that |o"i| > o"2 at intermediate heights, where |/3e| ^ 1 while 
/3i < 1. At higher z, typically |/3e| > A > 1 and the fluid is 
in the ambipolar diffusion regime. In these conditions, equa- 
tions (9) and (10) show that ai <^ (T2. On the contrary, near 
the rriidplane either component may dominate, depending 
on the fluid density and the strength of the magnetic field, 
ui will be larger if the fluid density is relatively small (larger 
radius) and/or the field is stronger. Conversely, if the field 
is weak and the density is sufflciently high, CT|| --^ (72 ^ cti 
and the conductivity will be resistive. This behaviour of the 
conductivity components as a function of height explains, in 
turn, the dependency of Xcrit with z: For weak fields it is 
<^ 1 near the midplane {<72 ^ <ti), but if B is sufficiently 
strong, Xcrit close to the midplane ^ 1. It then increases 
with z and levels in the intermediate sections of the disc at 
~ 1 (cTi » cr2)- Finally, Xcrit drops sharply near the surface 
(1T2 » fi again). 

We now briefly turn our attention to the dependency 
of the magnetic coupling with height. Note that x first in- 
creases with z, as expected, in response to the enhanced 



ionisation fraction, and lower fluid density, away from the 
midplane. Closer to the surface, it decreases again, as a re- 
sult of the abrupt drop in the magnitudes of the conductivity 
components. In the following subsections we explore more 
fully the dependency of x versus Xcrit with height, which 
indicates which non-ideal MHD regime is locally dominant, 
for the disc models of interest. 



4.3.1 Minimum-mass solar nebula disc 

Figs. 4 to 6 present curves of x (solid lines) and Xcnt = 
|o"i|/o"x (dashed lines) as a function of height for the fiducial 
model. Results are shown for Ji = 1, 5 and 10 AU and 
different choices of the magnetic field strength. Cosmic rays 
are either included (top panel of each figure) or excluded 
from the disc (bottom panels). Bottom (solid) and leftmost 
(dashed) curves correspond to B = ImG in all cases. B 
changes by a factor of 10 between curves, except that for 
5 and 10 AU the top (and rightmost) curves correspond to 
the maximum field strength for which perturbations grow. 

We note that in all plots, the magnetic coupling near 
the midplane increases with B. In these regions of the disc, 
both Pi and /?e are typically <C 1 (except possibly when B 
is very strong) and from equation (11), 



~ eerie ,a r, N 



(66) 



which is independent of B. As a result, x ^ B^- On the 
contrary, near the surface the /3j can be very large and 



eerie {0i — Be) 1 



(67) 



Because of this, near the surface x is quite insensitive to 
changes in B, as evidenced in Figs. 4 to 6. 

Note that when cosmic rays are assumed to be excluded 
from the disc, there is a discontinuity in the curve of x vs. 
z/H at R — 1 and 5 AU. This discontinuity is caused by the 
drop in the ionisation fraction at the height where x-rays 
are not able to penetrate any further within the disc (see 
bottom panels of Figs. 4 and 5). It is, therefore, not present 
at 10 AU, because x-rays reach the midplane at this radius 
(see bottom panel of Fig. 6) . 

In general, when B increases, the region around the 
midplane where Hall diffusion is locally dominant (x < 
\ai\/a±, SW03) is reduced, as a result of the stronger mag- 
netic coupling. At 1 AU, when cosmic ray ionisation is in- 
cluded, this criterion is satisfied in the inner sections of the 
disc for all magnetic field strengths of interest (Fig. 4, top 
panel). Above this region, there is a relatively small section 
where ambipolar diffusion is locally dominant while near the 
surface, ideal-MHD holds for all magnetic field strengths. 

For R = 5 AU (Fig. 5, top panel). Hall diffusion dom- 
inates for z/H < 2.2 when B < 10 mC, but for 10 mC 
< S < 100 mG, ambipolar diffusion is locally dominant 
in the inner sections of the disc (|cri|/(Tx < x < 10). For 
stronger fields, the magnetic coupling is such that x > 10 
even at the midplane and the fluid is in ideal-MHD condi- 
tions over the entire cross-section of the disc. At this radius, 
the fluid is in the resistive regime only for very weak fields 
(B < 1 mG) and very high densities (uh > 10^^ cm~^; War- 
die, in prep.). Similarly, results at i? = 10 AU (Fig. 6, top 
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panel), show that ambipolar diffusion dominates in the in- 
ner sections of the disc for B < 10 rriG. For stronger fields, 
ideal-MHD is a good aproximation at all z. 

When cosmic rays are assumed to be excluded from the 
disc, the previous results at 1 AU (bottom panel of Fig. 4) 
are largely unchanged, except for those obtained with B = 1 
G, where Ifall diffusion dominates now for z/H ^ 1.8, up 
from ~ 1 in the previous CEise. This is due to the sharp fall in 
the magnetic coupling in the region which x-rays are unable 
to reach. At 5 AU, the Hall regime is now dominant near 
the midplane for all studied magnetic field strengths, while 
for R = 10 AU there is now a Hall dominated central region 
for B < 10 mG. It extends to z/H 1.3. When B < 100 
mG, ambipolar diffusion dominates near the midplane, but 
for stronger fields ideal MHD holds for all z. 

Finally, we observe that (as expected), the magnetic 
coupling at the midplane increases with radius in response 
to the higher ionisation fraction in the central sections of the 
disc at larger radii. In fact, Xo increases by 3 - 4 orders of 
magnitude between R — 1 and 10 AU. On the contrary, the 
coupling at the surface does not change as much with radius, 
decreasing only from ~ 22 to ~ 2.5 between the same radii. 



4-3.2 A more massive disc 

Fig. 7 displays curves of x and Xcrit as a function of height 
for a more massive disc, as detailed in section 2.2. The ioni- 
sation balance is calculated assuming that cosmic rays pen- 
etrate the disc and results axe shown for R = 1, 5 and 10 
AU for the same range of magnetic field strengths explored 
in the minimum-mass solar nebula model. 

Increasing the surface density causes the ionisation frac- 
tion near the midplane to drop sharply. As a result, the mag- 
netic coupling is drastically reduced in these regions at all 
radii. This is especially noticeable at 1 AU, where in this 
case x-rays are completely attenuated for z/H < 2.6. This, 
together with a very weak ionisation rate {Ccr 10~^** for 
Z/H « 2.6 and it is negligible at « = 0), causes x to be very 
low -and fairly independent of z- in this section of the disc. 

The weaker coupling at low z results in Hall diffusion 
being locally dominant over a larger cross-section of the disc 
in this model. We find that for R = 1 AU, there is again a 
Hall dominated region about the midplane for all studied 
B. This section now extends to z/H ~ 2.3 for B = 1 G, up 
from z/B ~ 1 in the minimum-mass solar nebula disc. For 
weaker fields, these regions are also larger. 

For R = 5 AU, Hall diffusion is now dominant near the 
midplane for B < 1 G. For comparison, in the previous case 
Hall conductivity was locally important for B < 10 mG. 
Furthermore, ambipolar diffusion dominates over a small 
cross-section of the disc for all but the strongest (B > 1 
G) magnetic field strengths at this radius. This contrasts 
with the previous case, where ideal MHD held at all z when 
B > 100 mG. Even at 10 AU, Hall diffusion now dominates 
near the midplane when B < 10 mG. When B is stronger 
than this, but still < 100 mG, ambipolar diffusion predomi- 
nates for z/H < 2, while in the minimum-mass solar nebula, 
ideal-MHD held for B > 10 mG. 

In the next section we analyse the effect of different con- 
ductivity regimes being dominant at different heights, in the 
structure and growth rate of MRI perturbations. Different 
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Figure 3. Example of the dependency of the conductivity com- 
ponents parallel and perpendicular (j\ and a^) to the field, 
magnetic coupling x and Xcrit = l^ll/ci as a function of height. 
Results shown correspond to the fiducial model for R = \ AU 
and B = 10 mG. Note that for z/i/ < 2, <T|| ~ iT2 3> ci and 
the conductivity is resistive. There is then a central section for 
which |<Ti| 3> <T2 (the fluid is in the Hall conductivity regime), 
while for higher vertical locations ~S> ox and ambipolar 

diffusion dominates. In the region where x < Xcrit , Hall diffusion 
still dominates the structure and growth rate of the MRI (SW03). 

disc models, radial locations and magnetic field strengths 
are discussed. 



5 MAGNETOROTATIONAL INSTABILITY 
5.1 Structure of the perturbations 

5.1.1 All unstable modes at 1 AU 

Fig. 8 shows the structure and growth rates of all unsta- 
ble MRI perturbations for the fiducial model at 1 AU and 
B = 100 mG. Solid lines denote 5 Br while dashed lines cor- 
respond to SB^. This notation is observed in this paper in 
all plots that display the structure of the instability. 

The fastest growing perturbation in this case grows at 
f = 0.5020 and there arc 15 unstable modes with 0.1151 < 
v < 0.5020. Slow growing perturbations, with v < 0.1806, 
are active even at the midplane, while faster modes exhibit 
a central magnetically inactive (dead) zone (Gammie 1996, 
Wardle 1997). Moreover, the extent of the dead zone in- 
creases with the growth rate and for the fastest growing 
mode it extends to z/H ~ 1.6. We observe that 5Br and 
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Figure 4. Comparison of the local magnetic coupling x (solid 
lines) and Xcrit = |fi|/f± (dashed lines) for the fiducial model 
at -R = 1 AU and different choices of the magnetic field strength. 
In each case, Hall diffusion is dominant in the regions where x < 
Xcrit (SW03). Ambipolar diffusion dominates where Xcrit < X J; 
10. When x stronger than this, the fluid is in nearly ideal- 
MHD conditions (W99). From top to bottom (solid lines), and 
right to left (dashed lines), the magnetic field drops from 1 G to 
1 mG. B changes by a factor of 10 between curves in all cases. 
Top panel: Cosmic rays are present. Hall diffusion dominates for 
z/H < 3.5 to z/H < 1 for B increasing from 1 mG to 1 G, 
respectively. Bottom panel: Cosmic rays are excluded from the 
disc by protostellar winds. Note the discontinuity in x ^■t the 
height below which x-rays are completely attenuated. 



SBif,, specially in the slow growing modes, are asymmetrical 
about zero and their averages over vertical sections of the 
disc have opposite signs. This appears to be related to the 
dependency of these two fluid variables in equation (46). 

5.1.2 Most unstable modes at different radii, including 
cosmic ray ionisation 

The weakest strength for which unstable modes could bo 
computed was 1 mG. This is a computational, rather than a 
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Figure 5. As per Fig. 4 for il = 5 AU. Top (and rightmost) 
lines correspond here to B = 795 mG (top panel) and S = 615 
mG (bottom panel), which are the maximum field strengths for 
which perturbations grow in each scenario. When cosmic rays are 
included, the Hall regime dominates near the midplane for _B < 10 
mG. Similarly, for B < 100 mG, there is an intermediate region 
where ax/a^ < x < 10: ^^d ambipolar diffusion is dominant. 
For stronger fields, ideal-MHD holds for all z. When cosmic rays 
are excluded, Hall diffusion dominates near the midplane, and 
ambipolar diffusion is dominant at intermediate heights, for all B 
of interest here. 



physical limit: the dependency of Vmax with B in the weak- 
field limit, shows no evidence for a minimum required field 
strength for the instability to grow (see section 5.2.1). A 
comparison of the structure of the most unstable modes for 
the fiducial model at -R = 1, 5 and 10 AU is shown in Fig. 9. 
We display solutions for B from 1 mG, up to the maximum 
field strength for which unstable modes grow for each radius. 
Unless stated otherwise, this criterion is followed in all plots 
that show the structure of the fastest growing modes as a 
function of B in this study. Note that the maximum value 
of B is dependent on the disc model, the ionising agents 
incorporated and radius of interest. This does not imply 
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Figure 8. Structure of all unstable modes for the fiducial model at i? = 1 AU and B = 100 mG. In each case solid lines show 5Br while 
dashed ones correspond to 5B^ . The growth rate is indicated in the lower right corner of each panel. There are 15 unstable modes, with 
0.1151 < V < 0.5020. Some perturbations grow at the midplane, particularly when v < 0.1806. Modes growing faster than this show a 
dead region around the midplane, whose vertical extent increases with the growth rate of the perturbations. Note the asymmetry of SBr 
and SBj, about zero, especially in slowly growing modes. 



that 1 mG is necessarily the weakest field strength for which 
unstable modes exist. In fact, from 

Both the structure and growth rate of these perturba- 
tions are shaped by the competing action of different conduc- 
tivity components, whose relative importance change with 
height. We defer the analysis of their growth rate to sec- 
tion 5.2. Here, we discuss the structure as a function of z. 
Note first of all that at each radius, the wavelength of the 
perturbations increase with the magnetic field strength, as 
expected by both ideal-MHD and non-ideal MHD local anal- 
yses (Balbus & Hawley 1991, W99). 

At 1 AU (leftmost column of Fig. 9), for B < 500 mG, 



the region next to the midplane is a magnetically inactive 
zone (Gammie 1996, Wardle 1997). The extent of this region 
decreases as the field gets stronger. For example, for B — 1 
mG it extends from the midplane to z/H ~ 1.8, but when 
B > 500 mG, there is no appreciable dead region, a result 
of the stronger magnetic coupling close to the midplane and 
the relatively small wavenumber of the perturbations. At 
this radius, when B is relatively strong {B > 500 mG), the 
amplitude of these modes increase with z, a property that 
is typical of MRI perturbations driven by ambipolar diffu- 
sion. Ambipolar diffusion modes have this property because, 
as the local analysis (W99) indicates, in this limit the local 
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Figure 6. As per Figs. 4 and 5 for ii = 10 AU. Prom top to bot- 
tom (and right to left), lines correspond here to B = 100, 10 and 
1 mG. In this case, if cosmic rays are present, Hall diffusion is not 
dominant locally for any field strength, given the strong magnetic 
coupling near the midplane. Ambipolar diffusion dominates when 
_B < 10 mG, but for stronger fields idcal-MHD holds over the 
entire cross-section of the disc. If cosmic rays are removed, the 
Hall regime is relevant near the midplane only for B < 10 mG. 
When B > 100 mG, ideal MHD holds for all z. ~ 



growth of unstable modes increases with the magnetic cou- 
pling, which (except in the surface regions) increases with z 
(sec Fig. 3). As a result, the local growth rate of the MRI 
also increases with height in this regime and is able to drive 
the amplitude of global unstable modes to increase. This 
explains the shape of the envelope of these modes. Finally, 
note that when the magnetic field is wealc, the perturba- 
tions' wavcmimbcr is very high, and SBr and 5B^ are not 
symmetrical about z — 0. 

At 5 AU, when B > 100 mG, perturbations grow even 
at the midplane. For weaker fields, they exhibit only a very 
small dead zone which extends to z/H < 0.5. Note also that 
the envelope of these modes, particularly for B > 100 mG, 
is fairly flat. This is explained by recalling that the magnetic 
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Figure 7. As per Figs. 4 to 6 for R = 1, 5 and 10 AU for a 
more massive disc, such that = lOSo and = lOpo- Cosmic 
rays are present. From top to bottom (and right to left), lines 
correspond to B = 1, 0.1, 10^^ ^nd 10"'^ G {R = 1 and 5 AU), 
and B = 0.1, 10"^ and 10-3 G (i? = 10 AU). The increased 
surface density causes the magnetic coupling to drop near the 
midplane. As a result, Hall conductivity dominates over a larger 
section of the disc around the midplane, while ambipolar diffusion 
is important (at higher z) for stronger fields than it was the case 
in the minimum-mass solar nebula model. 
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coupling at the midplane in this region is very high (x ~ 90 
for B — 500 mG, Fig. 5, top panel), so the ideal-MHD ap- 
proximation holds. Under these conditions, unstable modes 
peak at the node closest to the midplane, given that the local 
growth rate is not a strong function of x and does not vary 
by much with height (see also top row of Fig. 6 of SW03, 
which shows similar envelopes for perturbations obtained 
with Xo = 100 and different configurations of the conductiv- 
ity tensor). On the other hand, for B < 100 mG, non- ideal 
MHD effects are important. When B = 10 mG, ambipolar 
diffusion is dominant for z/H < 2.3. This conductivity term 
is likely to drive this perturbation's structure, as evidenced 
by the central dead zone and the envelope peaking at an in- 
termediate height. For B = 1 mG, Hall conductivity is dom- 
inant close to the midplane {z/H < 2.2), the region where 
the envelope of this perturbation peaks. Furthermore, the 
high wavenumber suggests that its structure is determined 
by local effects. As a result. Hall conductivity is likely to 
drive the structure of this mode. 

Finally, for _R = 10 AU there is no appreciable dead 
zone for B > 10 mG, given the strong magnetic coupling 
at this radius. Even for B = 1 mG, Xo ~ 3, a figure that 
increases to ~ 100 for B = 100 mG. As a result, the fluid 
is in nearly ideal-MHD conditions (see section 4.3.1), which 
explains the flat envelope of these modes. 

5.1.3 Most unstable modes at different radii, excluding 
cosmic ray ionisation 

We also calculated the structure of the most unstable modes 
at the same three radial positions under the assumption 
that cosmic rays are excluded from the disc by protostellar 
winds (Fig. 10). In this case the central dead zones at 1 AU 
{B < 100 mG) extend over a larger cross-section of the disc, 
given that cosmic rays (when present) are the main source 
of ionisation near the midplane at this radius (see section 
4.2). Without them, the electron fraction below z/H ~ 1.7 
(where x-rays are completely attenuated) plummets, caus- 
ing the amplitude of MRI perturbations in this section of 
the disc to be severely reduced. As a result, when B < 100 
mG, perturbations are damped for z/H < 2. Another differ- 
ence with the previous case is noticeable for stronger fields 
(B > 500 mG), where before the MRI was active even at 
the midplane. In this case, the abrupt change in the ionisa- 
tion balance at the height where x-ray ionisation becomes 
active, (Fig. 1, top panel), causes the current to be effec- 
tively discontinuous there and produces the observed kink 
in the amplitude of these modes. 

In the zone where x-rays are excluded {z/H < 1.7), 
the main ionising agent is the decay of radioactive elements 
within the disc. However, they can only produce very weak 
magnetic coupling. For example, Xo ~ 10^* for B = 1 mG, 
and it only increases to ~ 10~^ when _B = 1 G. As a result, 
the amplitude of the perturbations near the midplane, even 
for strong fields {B > 500 mG), is very small at this ra- 
dius. If the abundance of metals in the gas phase is reduced 
from the fiducial value adopted here {S2 ~ 0.02) to (say) 
~ 2 X 10"^, the zone ionised only by radioactivity becomes 
magnetically dead (as expected). The perturbations' growth 
is only marginally affected. 

For R = 5 AU, we also observe a small kink in the 
perturbations' amplitude for B > 100 mG. In this case it 



occurs much closer to the midplane, because at this radius 
x-rays penetrate to z/H ^ 0.3. Finally, for _R = 10 AU, x- 
rays penetrate the entire cross section of the disc, so there 
is no kink in the amplitude of these modes. However, the 
midplane cosmic ray ionisation rate at this radius is about 
two orders of magnitude larger than that of x-rays (in fact 
they dominate over x-rays for z/H < 2.5, see Fig. 1, bot- 
tom panel), so excluding them does reduce significantly the 
magnetic coupling in this region. As an illustration of this, 
note that for B = IQ mG, Xo decreases from ~ 16 in the 
previous case to only ~ 1 here (see Fig. 6). As a result, for 
z/H < 0.3 the fluid is in the Hall regime while ambipolar 
diffusion dominates at higher z. Being a high wavenumber 
perturbation, the structure of this mode reflects mainly local 
fluid conditions. This explains the shape of its envelope (see 
rightmost column of Fig. 10, second panel from the top): 
flat envelope close to the midplane where the Hall effect is 
dominant and amplitude increasing with z at higher vertical 
locations, driven by ambipolar diffusion. On the other hand, 
for B = 1 mG, X < kil/cr± for z/H < 1.3. This, together 
with the high wavenumber of the perturbations, causes Hall 
diffusion to shape the envelope of this mode. 

5.1.1 Conductivity regime comparison (aiBz > 0) 

It is interesting to compare the structure of the most un- 
stable perturbations obtained with different configurations 
of the conductivity tensor (assuming that different conduc- 
tivity regimes dominate over the entire cross section of the 
disc) against the full conductivity results discussed in sec- 
tion 5.1.2. This way we can explore more fully the effects of 
different conductivity components in the overall properties 
of the instability. Fig. 11 presents such solutions for the fidu- 
cial model at _R = 1 AU as a function of the strength of the 
magnetic field. The left column shows solutions computed 
with a full conductivity tensor, while the middle and right 
ones display modes obtained using the ambipolar diffusion 
(iTi = 0) and Hall regime {(T2 = 0, aiB^ > 0) approxima^ 
tions, respectively. 

Note that, for all magnetic field strengths shown here, 
perturbations computed using the ambipolar diffusion ap- 
proximation have a more extended central dead zone than 
modes incorporating a full conductivity tensor. This is in 
agreement with the local analysis ( W99) , which showed that 
the growth of MRI perturbations in the ambipolar diffusion 
limit decreases steadily with the local magnetic coupling. 
As a result, the local growth of these modes is severely re- 
stricted near the midplane. This is especially effective for 
weaker fields, when global effects are less important, due 
to the high wavenumber of the perturbations. For exam- 
ple, when B = 1 mG, the magnetically dead zone in modes 
found using the ambipolar diffusion approximation extends 
to z/H ^ 2. The thickness of the dead region decreases for 
stronger B, but even with B = 1 G, there is a small section 
{z/H < 0.5) in which perturbations do not grow. 

Turning now our attention to Hall perturbations, we 
observe that their amplitude is fairly stable, especially for 
strong fields {B > 500 mG), where they show the charac- 
teristic flat envelope reported by SW03. When the field is 
weaker than this, perturbations are damped near the mid- 
plane. This is particularly evident for B < 100 mG and 
is probably related to the fact that at 1 AU the magnetic 
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Figure 9. Comparison of the structure of the most unstable modes for the fiducial model at ii" = 1, 5 and 10 AU for different choices of 
the magnetic field strength. B is indicated in the upper right hand of each panel while the growth rate is shown in the bottom right hand 
side. Results are shown for B spanning from 1 mG, the weakest magnetic field strength for which unstable modes could be computed, to 
the maximum strength for which unstable modes grow, which changes with the radius. Note that for all three radii, the wavelength of 
the perturbations increases with the strength of the magnetic field, as expected. At 1 AU, for B < 500 mG, there is a central dead zone 
which extends to z/H < 2. At 5 AU, only modes for B < 100 mG exhibit a dead region. At this radius the magnetic coupling, even at 
the midplane, is sufficiently large for ideal MHD to be a good approximation for stronger fields. Finally, for R= 10 AU, ideal MHD holds 
throughout the cross-section of the disc for S > 10 mG and perturbations in this region of parameter space grow even at the midplane. 
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Figure 10. As per Fig. 9 assuming cosmic rays are excluded from the disc. Note the increased extent of the central dead zone at 1 AU 
{B < 100 mG) and the kink in the amplitude of the perturbations for stronger fields. This feature is caused by the sharp change in the 
ionisation balance of the disc at the height where x-ray ionisation becomes active. As expected, this effect is less pronounced at 5 AU 
and unnoticed at 10 AU, where x-rays reach the midplane. 
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coupling at the midplane in this limit is very low. Even for 
B — 100 mG, Xo in the Hall limit ^ 10^^ and it d(>crc;ascs 
to ~ 10^^'^ for B = 1 mG (note that these values are sig- 
nificantly smaller than the ones shown in Fig. 4, which were 
obtained with a full conductivity tensor) . Local results high- 
light that the maximum local growth rate of MRI pertur- 
bations in the Hall limit is effectively unchanged from the 
ideal case when x ~* (W99). However, the global analysis 
presented here shows that when the coupling is very low, the 
amplitude of global Hall limit perturbations can be damped 
close to the midplane (see rightmost column of Fig. 11, top 
three panels). On the other hand, there is only a small dead 
zone for i3 ~ 1 mG, which extends from the midplane to 
z/H ~ 0.5. For stronger fields, there is no appreciable dead 
region. 

In general, the structure of full cr perturbations reflect 
the contribution of Hall as well as ambipolar diffusion con- 
ductivity terms. When the magnetic field is strong (B > 100 
mG), ambipolar diffusion is locally dominant over a more 
extended section of the disc, and the dead zone of perturba- 
tions in this limit is much smaller than they are for weaker 
fields. On the other hand, modes in the Hall limit grow 
now even at the midplane and have a significantly higher 
wavenumber than ambipolar diffusion perturbations. This 
is reflected in the structure of full <T modes in this region 
of parameter space: Their envelope is shaped by ambipo- 
lar diffusion (the amplitude increases with height), but the 
wavenumber is higher and they grow closer to the midplane 
than pure ambipolar diffusion modes do. This reflects the 
contribution of the Hall effect. These results are also in 
agreement with similar trends found with illustrative cal- 
culations in SW03. 

It was discussed in section 4.3 that at this radius (1 AU), 
Hall conductivity is locally dominant near the midplane for 
all B of interest. The transition to the zone where ambipo- 
lax diffusion dominates occurs at a lower z for stronger fields 
(Fig. 4, top panel) . This analysis reveals that the Hall effect 
modifies the structure of MRI modes even when Hall diffu- 
sion is locally dominant only for a small section close to the 
midplane of the disc. In section 5.2.1 it will be discussed how 
the Hall effect also alters the growth of all fastest growing 
modes at this radius. 

At 5 AU, ambipolar diffusion is important in the inner 
sections of the disc for B < 100 mG (Fig. 5). It is expected 
that perturbations obtained using this approximation will 
be different from the corresponding full cr modes in this 
region of parameter space. Fig. 12 compares solutions ob- 
tained with different configurations of the conductivity ten- 
sor at this radius. Ambipolar diffusion modes (left column) 
are indeed different from full cr {aiBz > 0) ones (middle 
column) when B < 100 mG. For stronger fields, ideal-MHD 
holds, so the structure of MRI modes computed using dif- 
ferent configurations of the conductivity tensor are alike, as 
expected. 

5.1.5 Conductivity regime comparison (aiB^ < 0) 

So far we have discussed solutions obtained with a aiB^ > 
Hall conductivity term, which corresponds to the case where 
the magnetic field and angular velocity vectors of the disc 
are parallel ($1 ■ B > 0). Wo now explore how those results 
are modified when these vectors are antiparallel. 



As noted before, we explored in this case a reduced re- 
gion of parameter space. In particular, solutions at 1 AU 
could not be computed at all, as Xo ^ 2 at this radius, 
the limit below which all wavenumbers grow in this regime 
(W99, see also section 4.1). On the other hand, results at 10 
AU do not differ appreciably in this case from those obtained 
using a ajBz > conductivity, given the strong magnetic 
coupling throughout the cross-section of the disc at this ra- 
dial location. At 5 AU, however, Hall diffusion is important 
for relatively weak fields {B < 10 mG), and both sets of 
results should be different in this region of parameter space. 
Fig. 12 (middle and right columns) compares the structure 
of MRI unstable modes computed with a positive and neg- 
ative Hall conductivity at this radius, as a function of the 
magnetic field strength. When aiBz < 0, no results could 
be computed for B < 8 mG, as Xo < 2 in these cases. We 
find that, indeed, for B = 10 mG the most unstable mode 
computed with a negative Hall conductivity has a higher 
wavenumber, and a slower growth rate, than the correspond- 
ing mode with aiBz > 0. 

5.1.6 More massive disc 

To finalise the analysis of the structure of MRI perturba- 
tions, we describe now the properties of unstable modes in 
a more massive disc, as characterised in section 2.2. Results 
are displayed in Fig. 13. As expected, the larger column den- 
sity causes unstable modes to have a more extended central 
dead zone in relation to results in the minimum-mass solar 
nebula model. This is particularly noticeable for R = 1 AU, 
where there is now a magnetically inactive zone for B < 1 G. 
X-rays are excluded from the midplane at this radius (they 
can only penetrate up to z/H ^ 2.6 for this disc model) 
and the cosmic ray ionisation rate at 2: = is negligible. 
As a result, even for B = 1 G the magnetic coupling at the 
midplane is ~ 10"'' and perturbations are damped at low z. 
We observe here the same trends discussed in the analysis 
of the minimum-mass solar nebula model in relation to the 
structure of these modes. Ambipolar diffusion shapes the 
envelopes, particularly for strong fields, while Hall diffusion 
increases their wavenumber. Finally, note the kink on the 
amplitude of the B = 5.85 G mode, attributed to the sharp 
increase in the ionisation fraction at the height where x-ray 
ionisation becomes active. 

At 5 AU, there is also a more extended dead zone up 
to B ~ 100 mG, in comparison with that of perturbations 
obtained using the mininmm-mass solar nebula model. In 
that case, ideal-MHD conditions held throughout the cross- 
section of the disc for B = 100 mG, with Xo ~ 20 (see sec- 
tion 5.1.2). On the contrary, in the present model, Xo ~ 0.2 
for this magnetic field strength and ambipolar diffusion is 
locally dominant for 0.7 < z/H < 2 (see Fig. 7, middle 
panel). As a result, the amplitude of this perturbation in- 
creases with z, a typical behaviour of ambipolax diffusion 
driven MRI, instead of being flat as before. 

Finally, for i? = 10 AU, ideal-MHD conditions hold 
throughout the disc cross-section for B > 100 mG (Fig. 
7, bottom panel). Because of this, the envelopes arc flat in 
this region of parameter space. On the other hand, when 
B ~ 10 mG, Xo ~ 1-5 and ambipolar diffusion dominates up 
to z/H 2. As a result, this mode exhibits a small central 
dead zone and the envelope peaks at an intermediate height. 
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Figure 11. Comparison of the structure of tlic most unstable modes of the MRI for the fiducial model at i? = 1 AU as a function of the 
magnetic field strength for different configurations of the conductivity tensor. The field strength spans from 1 mG up to 5 G, which is 
the maximum B for which perturbations grow in the ambipolar diffusion regime. The left column shows solutions obtained with a full 
conductivity tensor. The middle and right columns corresponds to the ambipolar (cri = 0) and Hall ((72 = 0, criBz > 0) approximations, 
respectively. Note that when the magnetic field is weak (B < 100 mG), full conductivity perturbations have a higher wavenumber, and 
grow closer to the midplane, than modes in the ambipolar diffusion limit. This reflects the contribution of the Hall effect. For stronger 
fields, the shape of the envelope appears to be driven by ambipolar diffusion (amplitude increases with height), with the Hall term 
increasing the perturbations' wavenumber. 
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Figure 12. Comparison of the structure of the fastest growing modes for the fiducial model at i? = -5 AU, for different configurations 
of the conductivity tensor. The left column corresponds to the ambipolar diffusion approximation. The middle and right columns show 
solutions obtained with a full conductivity tensor when the magnetic field and angular velocity vectors of the disc are parallel (ui > 0) 
and antipajrallel {aiB^ < 0), respectively. Note that when B < 10 mG, full tr solutions obtained with Hall conductivity terms of opposite 
signs are different, as this component dominates near the midplane in this region of parameter space. On the other hand, ambipolar 
diffusion perturbations differ from full <t ones for D < 100 mG, given that for these magnetic field strengths <T2 dominates close to the 
midplane. For stronger fields, the fluid is in nearly ideal-MHD conditions and perturbations in all three configurations of the conductivity 
tensor are alike. 



When B ^ 1 mG, Xo ~ 0.05 and Hall diffusion is dominant 
for z/H < 2. The central dead zone is much reduced, and 
this perturbation grows closer to the midplane despite the 
weak magnetic coupling. 



5.2 Growth rate of the perturbations 

5.2.1 Fastest growing modes for different conductivity 
regimes 

The dependency of the growth rate of the most unstable 
modes with the strength of the magnetic field is shown in 
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Figure 13. Comparison of tlie structure of the most unstable modes of tlie MRI for a more massive disc, incorporating cosmic ray 
ionisation. We present results at i? = 1, 5 and 10 AU as a function of the strength of the magnetic field. The disc surface and mass 
density are = lOSo and = lOpo- For simplicity, it was assumed that the temperature, sound speed and scale height are unchanged 
from those of the minimum-mass solar nebula model. Note that these modes have a more extended central dead zone, especially for 
R = 1 AU, in relation to the fiducial model. The kink in the perturbation at -R = 1 AU for B = 5.85 G (leftmost column, bottom panel), 
is caused by the sudden change in the ionisation balance at the height where x-rays are not able to penetrate further within the disc. 
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Fig. 14 for different configurations of the conductivity ten- 
sor. Results correspond to the fiducial model at i? = 1 AU 
(top panel) and 10 AU (bottom panel). 

The overall dependency of the Vmax vs. B curve is typ- 
ically as follows: Vmax initially increases with B with a ~ 
power law dependency. It then levels at ~ 0.75, the majd- 
mum growth rate for idoal-MHD perturbations in kcplcrian 
discs. Eventually, the maximum growth rate decays sharply 
at a characteristic magnetic field strength at which the per- 
turbations' wavelength is ~ H, the scaleheight of the disc. 
Note that in the weak-field limit, the dependency of the i^max 
vs. B curve shows no evidence of a minimum field strength, 
below which modes do not grow. This characteristic shape 
of the fmax vs. B curve can be explained in terms of the 
dependency of the perturbations' growth rate with B when 
different conductivity components are dominant at different 
heights. 

Examining Fig. 14 (top panel), it is clear that for weak 
fields {B < 200 mG), perturbations computed with a full cr 
have a Umax VS. B dependency similar to that of solutions 
obtained using the ambipolar diffusion approximation. This 
indicates that this feature is driven by ambipolar diffusion. 
In this limit, the local maximum growth rate of MRI modes 
increases with the magnetic coupling (W99), which in turn, 
increases with B, except when l3i and /3e S> 1 (e.g. near the 
surface). As a result, Vmax should increase with B, as ob- 
served. On the contrary, when the magnetic field is stronger 
than ~ 200 mG, the growth rate of MRI modes obtained 
with a full conductivity tensor is practically identical to that 
of Hall limit modes. In this region of parameter space, the 
growth rate of global modes remains unchanged when the 
magnetic field gets stronger, as expected when Hall diffusion 
drives the growth of the instability. Eventually, B becomes 
strong enough that the fastest growing mode becomes ~ the 
scaleheight of the disc, and the growth rate rapidly declines. 
These results are in agreement with previous findings by 
W99 and SW03 . 

Note also that the maximum growth rate of perturba- 
tions in the Hall limit does not change appreciably with B, 
until they are damped for a sufficiently strong field, as ex- 
pected. In fact, Hall perturbations grow at ~ the ideal-MHD 
rate even for B ^ 1 mG. 

Hall diffusion is locally dominant at 1 AU in the lower 
sections of the disc for all magnetic field strengths for which 
unstable modes exist (Fig. 4, top panel). As a result, pertur- 
bations obtained with a full conductivity tensor grow signifi- 
cantly faster than modes found using the ambipolar diffusion 
approximation for all B. 

For R = 5 AU, the Hall effect increases the growth rate 
of full cr perturbations when B < 10 mG, the region of pa- 
rameter space where Hall diffusion is locally dominant near 
the midplarie (compare the growth rates in left and middle 
columns of fig. 12. Moreover, full cr modes computed with 
aiBz > (middle column) and aiB^ < (right column) 
grow at different speeds and have adifferent structure, as 
expected in the region where Hall diffusion dominates. 

Finally, Fig. 14 (bottom panel) displays the growth rate 
of the most unstable modes at 10 AU. Results were com- 
puted with a full conductivity tensor -but with opposite 
signs of fl • B-, as well as using the ambipolar diffusion 
approximation. In this case. Hall diffusion is locally unim- 
portant for all B studied here, so the sign of aiB^ does not 



affect the growth of the perturbations. On the other hand, 
at this radius, ambipolar diffusion is important for B < 10 
mG. This slows the growth of perturbations in the ambipo- 
lar diffusion limit when the magnetic field is weaker than 
this value. 



5.2.2 Fastest growing modes at diferent radii 

Finally, we compare the growth rate of the most unstable 
modes at different radii as a function of the strength of 
the magnetic field (Fig. 15). Three sets of results are dis- 
played. The first two correspond to the minimum-mass so- 
lar nebula disc assuming cosmic rays either penetrate the 
disc (top panel) or are excluded from it (middle panel). The 
last set presents results for the more massive disc model 
(bottom panel), incorporating cosmic ray ionisation. In all 
three cases, MRI perturbations grow over a wide range of 
magnetic field strengths. 

Note that the maximum field strength for which un- 
stable modes exist (within each panel), is weaker at larger 
radii. In ideal-MHD conditions, as well as when either the 
ambipolar diffusion or Hall (ctiBz > 0) conductivity regimes 
dominate, unstable MRI modes are damped when va/cs ~ 1 
(Balbus & Hawley 1991), which corresponds to A ~ J?. In 
the Hall {aiB^ < 0) limit, unstable modes have been found 
for va/cs ~ 3 (SW03). In any case, as both the gas density 
(eq. 26) and the sound speed (etj. 17) decrease with radius, 
the ratio of the Alfven to sound speed associated with a 
particular magnetic field strength increases with R, and as 
a result, the perturbations are damped at a weaker field for 
larger radii. 

The maximum magnetic field strengths for which MRI 
unstable modes grow, as well as the range for which Vmax ~ 
the ideal-MHD rate, are summarised in Table 1 for the radii 
and disc structures of interest here. Note that for the fiducial 
model, we obtained unstable modes at 1 AU for B < 8 G. 
The growth rate is of the order of the ideal-MHD rate {f = 
0.75) for 200 mG ^ B < 5 G. When cosmic rays are assumed 
to be excluded from the disc (middle panel), unstable modes 
are obtained at 1 AU only for B < 2 G. The corresponding 
range at 5 AU is only slightly reduced and it is essentially 
the same as before at 10 AU. This is consistent with our 
expectation that cosmic rays are a particularly important 
source of ionisation at 1 AU, where x-rays are excluded from 
the midplane. In this case, Vmax ~ 0.75 at 1 AU for 100 mG 
^B<1G. 

In the massive disc model (bottom panel of Fig. 15), 
the MRI is active at 1 AU for B < 6 G and u^ax ~ 0.75 for 
200 mG < B < 500 mG. For i? = 5 and 10 AU, MRI un- 
stable modes exist in this disc for stronger fields than in the 
minimum-mass solar nebula model. This can be explained 
by recalling that the larger mass and column density trans- 
lates into a larger gas pressure and a stronger equipartition 
magnetic field strength at the midplane. MRI modes are, 
as a result, damped for a stronger field than they are in 
the minimum-mass solar nebula model. On the contrary, for 
R = 1 AU, the range of magnetic field strengths over which 
perturbations grow does not change appreciably (see Table 
!)■ 
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Table 1. Magnetic field strengths for which MRI perturbations grow at 1, 5 and 10 AU for different disc models. For the minimum-mass 
solar nebula disc, two sets of results are shown, where cosmic rays are either present or excluded from the disc. For the massive disc 
model, only the former case is presented. Columns labeled 'Bmax' list the maximum value of B for which each disc supports unstable 
MRI modes. Similarly, the ~ 0.75)' columns specify the subset of these for which perturbations grow at nearly the ideal-MHD rate 
(u ~ 0.75). 
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Radius (AU) 


Bmax 
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Bmax 
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0.80 
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0.62 
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2.. 36 


0.05 - 0.5 


10 


0.25 


0.002-0.05 


0.25 


0.005 - 0.05 


0.82 


0.01 - 0.5 



6 DISCUSSION 

In this paper wo have explored the vertical structure and 
linear growth of MRI perturbations of an initially vertical 
magnetic field, using a realistic ionisation profile and assum- 
ing that ions and electrons are the sole charge carriers. This 
formulation is appropriate to model low conductivity discs 
(Gammie & Menou 1998; Menou 2000; Stone et al. 2000) at 
late stages of accretion, after dust grains have settled into a 
thin layer around the midplane (~ 10^ years, Nakagawa et. 
al. 1981, DuUemond & Dominik 2004) and become dynam- 
ically uncoupled from the gas at higher z. Solutions were 
obtained at 1, 5 and 10 AU from the central object as a 
function of the strength of the magnetic field for diS'erent 
configurations of the conductivity tensor, disc model and 
sources of ionisation. 

We have shown that the magnetic field is dynamically 
important in low conductivity accretion discs over a wide 
range of field strengths. An example of this activity is the 
generation and sustaining of MRI-unstable modes, which are 
needed to provide the angular momentum transport for ac- 
cretion to proceed. The structure and growth rate of MRI 
perturbations axe a function of the disc properties and the 
strength of the field. For a particular radius and disc model, 
they are a result of the competing action of different domi- 
nant conductivity components at different heights. 

For R ^ & few AU and relatively weak fields {B < 100 
mG at 1 AU or S < 1 mG at 5 AU), the fluid close to the 
midplane is in the Ohmic conductivity regime ((T|| ^ a2 
|c7"i|). Conversely, for larger radii (or stronger fields) Hall 
conductivity dominates even at 2 = (Wardle, in prepa- 
ration; see also section 4.3). Ambipolar diffusion dominates 
above this region until x 10 (W99). For still higher z, the 
magnetic coupling is so strong that ideal-MHD holds. The 
heights at which these transitions take place, for a particular 
disc model and radial position, are a function of the strength 
of the magnetic field. Note that even when \ai\ <§; (72, Hall 
diffusion still dominates the structure and growth of MRI 
modes when x Xcnt = \ai\/a± (SW03). 

At 1 AU, for the fiducial model, unstable modes exist 
for S < 8 G. For a significant subset of these strengths (200 
mG ^ B < 5 G), perturbations grow at about the ideal- 
MHD rate {— 0.75il in keplerian discs). At this radius. Hall 
diffusion is locally dominant near the midplane for all mag- 
netic field strengths for which MRI unstable modes exist. 
As a result, modes computed with a full conductivity tensor 
have a higher wavenumbcr, and grow faster, than perturba- 
tions obtained using the ambipolar diffusion approximation. 



On the other hand, ambipolar diffusion shapes the envelopes 
of these modes (especially for strong fields) , causing them to 
peak at an intermediate height, instead of having the char- 
acteristic flat envelope that modes in the Hall limit have for 
these field strengths. 

Cosmic rays are an important source of ionisation at 
1 AU, given that x-rays do not reach the midplane at this 
radius (see Fig. 1). Consequently, if they are assumed to be 
excluded from the disc, the extent of the magnetically inac- 
tive 'dead zone' about the midplane increases and unstable 
modes grow only for S < 2.1 G. Furthermore, under this 
assumption, perturbations obtained with a relatively strong 
initial magnetic field (e.g. B > 500 mG), exhibit a kink 
in their amplitude. This feature is attributed to the sharp 
increase in the ionisation fraction at the height where x- 
ray ionisation becomes active. Finally, for the massive disc 
model (Eq = lOSo and = lOpo, incorporating cosmic 
ray ionisation) , the central dead zone extends to a higher z 
at 1 AU, as the larger column density associated with this 
model causes the ionisation fraction at low z to drop sharply 
(the disc is shielded from x-rays below z/H ~ 2.6 and even 
the cosmic ray ionisation rate at the midplane is negligible). 
MRI perturbations grow in this case for B < 6 G, a range 
not significantly different from that of results in the fiducial 
model. This is the case because the effective surface den- 
sity is largely unchanged, given the shielding of the inner 
sections. 

At 5 AU, for the fiducial model, MRI modes grow for 
B < 800 mG. When 20 mG < B < 500 mC, the growth 
rate is ~ 0.75n. At this radius, Hall diffusion increases the 
growth rate and wavenumbcr of unstable modes for B < 100 
mG, but for stronger fields, perturbations obtained with dif- 
ferent configurations of the conductivity tensor have similar 
structures and growth rates, a signal that ideal-MHD condi- 
tions hold (Fig. 12). We conducted an analysis of the effect 
of the alignment of the magnetic field and the angular ve- 
locity vectors of the field. Results indicate that the sign of 
aiBz is important for B < 10 mG, the range of magnetic 
field strengths for which Hall diffusion is locally dominant 
at low z. As in the 1 AU case, excluding cosmic rays reduces 
the range of magnetic field strengths for which MRI modes 
grow (unstable modes are found for B < 600 mG in this 
case) and there is a kink in the amplitude of the perturba- 
tions for B > 100 mG. Finally, when the surface density of 
the disc is increased, the central dead zone occupies a larger 
cross-section, as expected, and unstable modes are found for 
stronger fields than in the fiducial model. This is a result of 
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Figure 14. Growth rates of the fastest growing modes as a 
function of the strength of the magnetic field for different con- 
figurations of the conductivity tensor. Results are presented at 
R = 1 AU (top panel) and 10 AU (bottom panel) for the fidu- 
cial model. It is evident that at 1 AU, when the magnetic field 
is weak (B < 200 mG), perturbations are driven by ambipolar 
diffusion. As a result, Umax increases with B. As the field gets 
even stronger, the maximum growth rate of full cr and Hall limit 
perturbations arc practically identical, signalling that Hall con- 
ductivity determines the growth of global unstable modes in this 
region of parameter space. Finally, Umax decreases rapidly when 
B is so strong that the wavelength of the fastest growing mode is 
~ the scaleheight of the disc. 



the stronger equipartition magnetic field strength for this 
disc model at this radial position. 

At 10 AU, the MRI is active for B < 250 mG and 
the growth rate is close to the ideal-MHD rate for 2 mG 
< B < 50 mG. Furthermore, for B < 10 mG, perturbations 
obtained with a full conductivity tensor grow significantly 
faster than modes in the ambipolar diffusion limit, which 
refiects the contribution of the Hall effect (Fig. 14, bottom 




-3 -2-10 1 
Log(B) 



Figure 15. Growth rate of the most unstable modes of the MRI 
for R = 1, 5 and 10 AU as a function of the strength of the 
magnetic field. Top and middle panels show results obtained with 
the fiducial model, assuming cosmic rays are either included, or 
excluded, from the disc. Bottom panel displays solutions obtained 
with the more massive disc, including cosmic ray ionisation. 
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panel). At this radius, when cosmic rays are excluded, the 
range of magnetic field strengths for which unstable modes 
exist is not affected, given that x-rays are able to penetrate 
to the midplane. In the more massive disc, the extent of 
the dead zone increases, especially for weak fields {B < 10 
mG) and perturbations grow for stronger fields than in the 
minimum-mass solar nebula disc. 

In the present formulation, it was assumed that ions and 
electrons are the only charge carriers. As discussed earlier, 
this is valid in late evolutionary stages of accretion, when 
dust grains have sedimented enough towards the midplane 
that they can be neglected when studying the dynamics of 
the gas at higher z. The timescale for this settling to occur 
is expected to be affected by MHD turbulence. Although in 
quiescent discs dust grains may quickly settle into a thin 
sub-layer about the midplane, the vertical stirring caused 
by MHD turbulence could potentially transport them back 
to higher vertical locations, preventing them from settling 
below a certain height (DuUemond & Dominik 2004 and ref- 
erences therein). The efficiency of this process is dependent 
on the disc being able to support MHD turbulence in the ver- 
tical locations where dust grains are present. On the other 
hand, dust grains can reduce the abundance of free elec- 
trons, and the efficiency of MHD turbulence itself, by pro- 
viding additional recombination pathways on their surfaces. 
To make things more complex, it is likely that the tran- 
sition between vertical sections where dust grains axe well 
mixed with the gas phase and those completely depleted of 
grains by settling, occur over a finite thickness (DuUemond 
& Dominik 2004). Ultimately, the equilibrium structure of 
accretion discs will reflect the complex interplay between all 
these processes. 

It is expected that the structure and growth of MRI 
unstable modes in such environments will be afi'ected by the 
dependency of the ionisation balance with height in the pres- 
ence of chemistry taking place on grain surfaces. The study 
of the properties of the MRI in a disc where dust dynamics 
and evolution are determined consistently, and where Hall 
conductivity is taken into account, is essential to understand 
more fully the presence and efficiency of MRI-driven angu- 
lar momentum transport in accretion discs. Nevertheless, 
from the results presented in this paper, it is clear that Hall 
diffusion is crucial for the realistic modelling of the magne- 
torotational instability in protostellar discs, particularly at a 
distance of ~ a few AU from the central object. More gener- 
ally. Hall conductivity is an important factor when studying 
the magnetic activity of low conductivity discs at these radii. 



7 SUMMARY 

We have presented in this paper the vertical structure and 
linear growth of the magnetorotational instability (MRI) in 
weakly ionised, stratified accretion discs, assuming an ini- 
tially vertical magnetic field. Both the density and the con- 
ductivity are a function of height and the conductivity is 
treated as a tensor and obtained with a realistic ionisation 
profile. Two disc models were explored: The minimum-mass 
solar mebula disc (Hayashi 1981, Nakagawa & Nakazawa 
1985) and a more massive disc, with the mass and surface 
density increased by a factor of 10. This formulation is ap- 
propriate for the study of weakly ionised astrophysical discs. 



where the ideal- MHD approximation breaks down (Gammie 
& Menou 1998; Menou 2000; Stone ct al. 2000). The ionisa- 
tion sources relevant here, outside the inner 0.1 AU from the 
central object, are non-thermal: Cosmic rays, x-rays and ra- 
dioactive decay. For the minimum-mass solar nebula model 
we compare solutions obtained including all three sources of 
ionisation with those arrived at assuming that cosmic rays 
are excluded from the disc by the protostar's winds. Recom- 
bination processes are taken to occur in the gas-phase only, 
which is consistent with the assumption that dust grains 
have settled into a thin layer about the midplane, and ions 
and electrons are the only charge carriers. Perturbations of 
interest have vertical wavevectors (fe = kz) only, which are 
the most unstable modes (when initiated from a vertically 
aligned magnetic field) in both the Hall and Ohmic regimes 
(Balbus & Hawley 1991; Sano & Miyama 1999). This is not 
necessarily the case in the ambipolar diffusion limit (Kunz 
& Balbus 2003), where the fastest growing modes can have 
radial as well as vertical wavenumbers. Under the adopted 
approximations, the properties of the MRI in the ambipolar 
diffusion and Ohmic conductivity limits are identical. 

Three parameters were found to control the dynamics 
and evolution of the fluid: (i) the local ratio of the Alfven to 
sound speed (va/cs); (ii) The local coupling between ionised 
and neutral components of the fluid (x), which relates the 
frequency at which non-ideal effects are important with the 
dynamical (keplerian) frequency of the disc; and {iii) the 
ratio of the components of the conductivity tensor perpen- 
dicular to the magnetic field (ai/cr2), which characterises 
the conductivity regime of the fluid. These parameters were 
evaluated at i? = 1, 5 and 10 AU for a range of magnetic 
field strengths. The linearised system of ODE was integrated 
from the midplane to the surface of the disc under appro- 
priate boundary conditions and solutions were obtained for 
representative radial locations of the disc as a function of 
the magnetic field strength and for different configurations 
of the conductivity tensor. 

The main results of this study are summarised below. 

(i) For the minimum-mass solar nebula model, incorpo- 
rating cosmic ray ionisation (the fiducial model): 

At 1 AU, unstable MRI modes exist for i? < 8 G. When 

200 mG ^ B < 5 G, the most unstable modes grow at ~ 
the ideal- MHD rate (= 0.75fi). Hall diffusion dominates the 
structure and growth rate of unstable modes for all magnetic 
field strengths for which they grow. For strong fields, am- 
bipolar diffusion shapes the envelope of the perturbations, 
which peak at an intermediate height. Finally, at this radius, 
a magnetically dead zone (Gammie 1996, Waxdle 1997) ex- 
ists when B < 1 G. As expected, the vertical extent of this 
zone decreases when the magnetic field gets stronger and/or 
when a more massive disc model is used. 

At 5 AU, MRI modes grow for B < 800 mG and the 
growth rate is close to the ideal MHD rate for 20 mG < B < 
500 mG. Perturbations incorporating Hall conductivity have 
a higher wavermmbcr and grow faster than solutions in the 
ambipolar diffusion limit for B < 100 mG. Unstable modes 
grow even at the midplane for B > 100 mG but for weaker 
fields, a small dead region exists. 

At 10 AU, the MRI is active for B < 250 mG. The growth 
rate is close to the ideal-MHD rate for 2 mG ^ B < 50 mG 
and when B < 10 mG, perturbations obtained with a full 
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conductivity tensor grow significantly faster than modes in 
the ambipolar difFusion hmit. Modes show only a very small 
dead region when B ~ 1 mG. 

(ii) When the magnetic field is weak (e.g. B < 200 mG at 
1 AU), the maximum growth rate of unstable MRI modes 
(i^max) increases with the strength of the magnetic field, a 
feature driven by ambipolar difTusion. 

(iii) When cosmic rays are assumed to be excluded from 
the disc, unstable modes at 1 AU grow only for B < 2.1 G. 
Results at 5 AU only change slightly, while solutions at 10 
AU are not affected at all, as expected. 

(iv) For the massive disc model (E^, = lOSo and = 
Wpo, incorporating cosmic ray ionisation), MRI perturba- 
tions grow for stronger fields at 5 and 10 AU, in relation to 
the minimum-mass solar nebula model. Results at 1 AU are 
unchanged, as in this case, the effective surface density is 
not significantly different. 

MRI perturbations grow in protostellar discs for a wide 
range of fluid conditions and magnetic field strengths. Hall 
diffusion largely determines the structure and growth rate 
of perturbations at radii of order of a few AU from the cen- 
tral protostar. This indicates that, despite the low magnetic 
coupling, the magnetic field is dynamically important in low 
conductivity astrophysical discs and will impact the dynam- 
ics and evolution of these discs. 
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